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Recently, Ishii et al.EP measured the magnetic suscep- 
tibility x for Cu 3 Cl 6 (H 2 0) 2 • 2H 8 C 4 S0 2 , which is con- 
sidered to be a quasi-one-dimensional material consist- 
ing of S = | trimer spin chains. The result indicates 
that x vanishes in the low temperature limit. They also 
measured the magnetization process for this material, 
and showed that there is a plateau of zero magnetization 
below the critical field H c ~ 3.9 T. From these experi- 
mental results, they concluded that the ground state is a 
singlet state with spin gap. The spin gap A is estimated 
as A/fce — 5.2 K from the value of H c . 

The proposed HamiltoniartP representing a spin chain 
in this material is given by 

H = Ji (Saj-i ■ S 3 j + S 3j ■ S 3j+ i) 
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where Sj is the S = 5 spin on site j. Three spins Saj—i, 
S 3 j and S$j+i form a trimer. The lattice structure is 
shown in Fig. 1. Three kinds of exchange constants Ji, 
J 2 and J3 are inferred to be positive and to satisfy the 
relation Ji > J 2 , J3 from the lattice parameters of the 
materially Hereafter, we use the unit of J\ = 1. 

The symmetric,case of J3 = Ji(= 1) has been studied 
by Takano et al£P and the system has been called the 
diamond chain. They almost exactly showed that there 
exist three phases in the parameter space; the ferrimag- 




Fig. 1. The lattice structure of the distorted diamond chain rep- 
resenting Cu3Clg(H20)2 ■ 2HgC4S02- Three kinds of exchange 
constants are shown. 




Fig. 2. Spin gaps of finite systems with L = 12, 18 and 24 (solid 
circles) for several values of J3 at J2 = 1. Each solid line repre- 
sents eq. (|^) with fitting parameters ci, C2 and Aoo for a value 
of J 3 . 



netic phase for J 2 < 0.909, the tetramer-dimer (TD) sin- 
glet phase for 0.909 < J 2 < 2 and the dimer-monomer 
(DM) singlet phase for J 2 > 2. The TD phase is a disor- 
dered phase with spin gap which originates from frustra- 
tion among exchange interactions, while the DM phase 
is a spin fluid phase without spin gap. Okamoto et all]) 
studied the general case of J3 ^ 1; i. e. the distorted 
diamond chain. The three phases develop in the J2-J3 
plane. They numerically determined the phase bound- 
aries. Also Tonegawa et al.Er numerically studied the 
magnetization process and showed plateaux for i and | 
of the saturation field. 

In this article, we estimate the values of the spin gap 
by the numerical diagonalization. Then we produce a 
contour map in the J 2 - J3 parameter space. The contour 
map represents an overall feature of the gapped phase 



of the S 



Heisenberg model on the distorted dia- 



mond chain. When further experimental information on 
Cu 3 Cl 6 (H 2 0) 2 •2H 8 C4S0 2 is given, the contour map will 
be useful to determine the values of the exchange con- 
stants for the real material. 

We first calculate the spin gap A l for finite chains with 
system size L. The spin gap Aoo in the thermodynamic 
limit is evaluated by extrapolation. We assume the size 
dependence of Al as 

A L = A 00 + r j- + ^ 2 (2) 

with constants c\ and c 2 . The numerical diagonalization 
has been done for L =12, 18 and 24 under the periodic 
boundary condition. We determine ci, c 2 and Aoo by 
fitting. In Fig. 2, we show Al as a function of L for 
several values of J 3 at J 2 = 1. For J 3 = 0, the estimated 
value of Aqo is about 0.002 and is close to zero; t 
nonzero value is interpreted as an extrapolation error 
For < J 3 < 0.4, the true value of the spin gap is very 
small or may be regarded as zero, since the estimated 
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Fig. 3. The spin gap Aoo in the thermodynamic limit (symbols) 
as a function of J3 for J2 = 0.7, 1.0 and 1.4. Each solid line 
represents D(Js) (eq. (g)) with fitting parameters a\ and 02 for 
a value of J 2 ■ 



Fig. 4. The contour map for the spin gap in the gapped phase of 
the J2-J3 plane. Bold solid lines rare the phase boundaries that 
Okamoto et al. have determined. a 



values are less than 0.002 and are within the extrapola- 
tion error. For J 3 > 0.6, the figure shows that the system 
has a finite spin gap. For J 3 = 0.5, Aoo is 0.0034, which 
is small but seems to be finite. This is consistent with 
the result of Okamoto et al. that the spin gap opens at 
the critical value Jf ~ 0.35 for J 2 = 1.eP 

In general, the spin gap in a dimer phase is expo- 
nentially small near the phase boundary to a spin fluid 
phase. Hence it is difficult to estimate A M near the 
boundary in the present case. To overcome this diffi- 
culty, we assume that J 3 dependence of the spin gap is 
given by 

D(J 3 )=a 1 ^jr~Jlexp(^ 1 ^ 7 ^j (3) 

for J 3 ~ J$ where a% and a 2 arc constants. We 
have the valuesjpf J| by inspecting the phase diagram of 
Okamoto et al.;H> e. g. J§ = 0.374, 0.354 and 0.460 for J 2 
= 0.7, 1.0 and 1.4, respectively. We carry out the fitting 
of the extrapolation data A M to eq. (|^) and determine 
ai and a 2 . Figure 3 represents the fitting function D(J 3 ) 
and the extrapolation data. We find that the extrapola- 
tion data are well reproduced by eq. (||) for A m > 0.02. 
Hence the function form in eq. ^ is reliable. We use 
eq. (g) to estimate the spin gap for A m < 0.02 near the 
critical value Jf . For example, the spin gap is estimated 
as 1.0 x 10" 2 at J 3 ~ 0.57, 1.0 x 10~ 3 at J 3 ~ 0.50, 
1.0 x 10~ 4 at J 3 ~ 0.47 and 1.0 x 10~ 5 at J 3 ~ 0.45 for 
■h = 1.0. 

Using these results, we draw contour lines of the spin 
gap in the J 2 -J 3 plain. The resultant contour map is 
shown in Fig. 4. We have calculated A m at the dis- 
crete positions (J 2 , J 3 ) with J 2 — 0.7, 0.8, 2.0 and 
J 3 = 0.5, 0.55, 1.0. For A m > 0.02, the positions of 
solid circles are determined by the linear interpolation 
among the spin gaps Aoo at the discrete positions. For 
Aqo < 0.02, the positions of open circles are determined 



by using -D(J 3 ) (eq. (||)) instead of Aoo. 

The temperature dependence of the experimental mag- 
netic susceptibility has a broad peak at ^70 K. It sug- 
gests that the energy scale of the characteristic exchange 
constant J\ is larger than 70 K. Here we consider a case 
of Ji being 100 K as an example-tP In this case, we have 
Aoo ~ 0.05 according to the observed spin gap ~5 K. 
Then J 2 and J 3 are limited to values close to the con- 
tour line of Aoo = 0.05 and of J 2 < 1. 
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